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Abstract 

We discuss neutrino oscillations in the framework of the quantum field theory with- 
out introducing the concept of neutrino weak eigenstates. The external particles are 
described by wave packets and the different mass eigenstate neutrinos propagate be- 
tween the production and detection interactions, which are macroscopically localized in 
space-time. The time-averaged cross section, which is the measurable quantity in the 
usual experimental setting, is calculated. It is shown that only in the extremely relativis- 
tic limit the usual quantum mechanical oscillation probability can be factored out of the 
cross section. 
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1 Introduction 



Neutrino oscillations have long been recognized as a powerful tool to probe the intrinsic 
properties of neutrinos 0]. Furthermore, it already appears that they may provide an 
elegant solution to the Solar Neutrino Problem, possibly leading to information on the 
basic properties of neutrinos such as mass and mixing angle. 

If neutrinos are massive and mixed, a weak eigenstate neutrino which is produced in a 
weak process accompanying a lepton is a linear superposition of mass eigenstates. In the 
standard treatment of neutrino oscillations [|l], 0] , the mass eigenstates are assumed to be 
relativistic and to have the same momentum but different energies. Because of the energy 
differences, the quantum mechanical probability of finding weak eigenstates becomes a 
function of the distance from the production point, leading to neutrino oscillations. 

Although the standard approach of treating neutrino oscillations with use of the weak 
eigenstates is physically intuitive and simple, it is, strictly speaking, neither rigorous nor 
sufficient for a complete understanding of the physics involved in neutrino oscillations. 
Furthermore, as shown in Ref. P], the usual "weak eigenstates" |z/q) = J2a ^aaWa) 
is the mixing matrix of the neutrino fields and {ua) are the mass eigenstates) describe 
correctly the neutrinos produced and detected in weak-interaction processes only in the 
extremely relativistic limit. Also, energy-momentum conservation in the process in which 
the neutrino is created implies that the different mass-eigenstate components must have 
different momenta as well as different energies [Q. On the other hand, if the particles 
involved in the production (as well as the detection) process are assumed to have def- 
inite four-momenta, then the neutrino produced (detected) is forced to have a definite 
four-momentum, implying that the neutrino is one of the mass eigenstates ||5[. This 
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observation suggests an apparent incompatibility between energy-momentum conserva- 
tion in the production and detection processes and the neutrino oscillations. In other 
words, energy- momentum conservation which forces neutrinos to be in mass eigenstates 
is incompatible with neutrino oscillations. This apparent incompatibility, however, does 
not arise in a physical situation since a necessary condition for neutrino oscillations to 
occur is that the neutrino source and detector are localized within a region much smaller 
than the oscillation length, and hence the neutrino momentum must have at least the 
corresponding spread given by the uncertainty principle 0. This spread is responsible 
for neutrino oscillations. 

The localization of the neutrino source and detector and the spread of the neutrino 
momentum imply that a propagating flavor neutrino is not described by a superposition 
of plane waves, but instead by a superposition of localized wave packets 0. The wave 
packet treatment necessary for a correct quantum mechanical description of neutrino 
oscillations has been discussed in Ref. |0. However, neutrino oscillations have, so far, been 
discussed in the framework of quantum mechanics of neutrino propagation, whereas the 
effects of the production and detection weak-interaction processes have been neglected. 
As shown in Ref. , the neutrino oscillation probability is independent from the details 
of the production and detection processes only in the extremely relativistic limit. Hence a 
quantum field theoretical treatment of neutrino oscillations is necessary for the discussion 
of the case in which some of the mass eigenstates happen to be not extremely relativistic. 
It is, of course, expected that a quantum field theoretical treatment must reproduce the 
quantum mechanical oscillation probability in the extremely relativistic limit. 

In this paper, we present the quantum field theoretical treatment of neutrino oscil- 
lations by using, as an example, a specific fiavor changing process (see Eq.(|l|)) and by 
calculating its cross section. The external (initial and final) particles are described by 
wave packets and the mass eigenstate neutrinos propagate from the production region 
to a detector which are macroscopically separated in space-time. Since energies and 
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momenta of the external particles are not precisely defined and energy-momentum is 
conserved in the interaction vertices within the uncertainty principle, the contributions 
from the propagation of different mass eigenstate neutrinos can interfere to produce os- 
cillations. 

In the quantum field theoretical treatment it is impossible to derive a general for- 
mula for neutrino oscillations because the cross sections depend on the details of the 
specific production and detection interactions involved. Hence in this paper we illustrate 
a general method using a specific example. Furthermore, strictly speaking, an oscillation 
probability cannot even be defined because the space dependence of the cross section can- 
not be factorized out. This does not imply that neutrino oscillation phenomena do not 
take place. The phenomena can only be inferred from actual measurements of cross sec- 
tions. It will be shown that the quantum field theoretical treatment yields the standard 
quantum mechanical oscillation probability only in the extremely relativistic limit. 

The plan of this paper is as follows: In Section |^ we present the calculation of the 
amplitude for a specific flavor changing process. In Section ^ we calculate the cross 
section and its time average, which is an experimentally measurable quantity. 

Finally, we discuss, in Section the case of extremely relativistic neutrinos and 
reproduce the standard oscillation probability. 



2 Amplitude 

Let us consider the weak flavor-changing processes 

^ ^u + Dj^ Dp + e- , 

occurring through the intermediate propagation of a neutrino, where Pj and Pp [Dj and 
Dp) are the initial and final production (detection) particles. For simplicity, we denote 
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the energy- momentum four vectors for the particles involved in Eq.([l|) as 



P = P' + PfM + Pu , Pu + k = k' + Pe 



(2) 



We will consider a process in which the production and detection interactions are macro- 
scopically localized at the coordinates (Xp,Tp) and (X£i,Td), respectively. 
The relevant weak interaction Lagrangians are 



G 



X] 



(3) 



where (x) and J^ix) are the weak currents of the production and detection particles, 
respectively and the other notations are self-evident. 
The amplitude for the process is 



T 



J d^xi d'^X2 Cp{xi)Cr){x2) 



(4) 



where the initial and final particles are described by the wave packets 



(27r)3/2 



^P,(p;Xp,Tp,(p-)) \Pj{p)) 



(5) 



i)e{Pe;XD,TD,{Pe)) 6 {p. 



(27r)3/2 

The form of the wavefunctions in momentum space is determined by how the initial 
particles are prepared and how the final particles are detected. In the following we 
will assume, for simplicity, gaussian wavefunctions, whose form is given by Eq.(p9D in 
Appendix A. In Eq.(^), (p) , ■ ■ ■ , (pe) are the average momenta of the particles around 
which their momenta are spread due to the uncertainty principle. The wave packets are 
constructed in such a way that at the time t = Tp the wave packets of the muon and the 



5 



production particles overlap at x ~ Xp and at the time t = Td the wave packets of the 
electron and the detection particles overlap at x ~ X^- 
The propagators of the mass eigenstate neutrinos are 



Hence the amplitude in Eq.(|p becomes, with Eqs. (]^) and 



(6) 



dk 



-iq-{x2-xi) 



X 



(27r)3/2 



„/, r^—ik x2 



(2vr)3/2 J (27r)3/2 



(7) 



X 



(27r)3/2 



(27r)3/2 
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where (A;, /c') and Jg{p,p') are the matrix elements of the weak currents of the detec- 
tion and production particles, respectively and we have not explicitly written down the 
arguments of the wave packets. 

The momentum integrations of external particles can easily be carried out if the wave 
packets in momentum space are sharply peaked around their average momenta (we have 
assumed, for simplicity, that all the wave packets are gaussian with the same width 
ax). After straightforward but tedious integrations, the amplitude can be written, with 
appropriate changes of the coordinates, as 

dS 



a 

X J d^xi exp 
X / d'^X2 exp 



(2vr) 



D 



q2 _ ^2 _|_ 



Vp exp —iqoT + iq-L 



-t {Ep -q,)t, + i{pp-q)-Xi-3^ + 3^ti- ^ tf 



-i [Ed + qQ)t2 + i [pd + q) ■ X2 - 3 — + 3 — — 12 

4o-H 2al 



W+u"-+v, 



where we have defined 



Pp= {P) - {P) - (Pm) 
3wp = V + + 

17^^ ji'((fc),(fc'))lZ;((Pe))7"(l + 75) 



Eo = {ko) - (K) - {E,) 
Pry ^ (k) - (k') - (p,) 

3V£) = U + u' + Ve 

Vp^{l+^,)rv,{{p,))Jn{p),{p')) ■ 



(9) 



Also in Eqs-d^) and 



we have introduced the notation 
(P) -/ (P') - 



u 



iPo) 

(k) 



u 



(P'o) 

{k') 



{Ee) 



iko) ' {k'o) 
Carrying out the integrals over Xi and X2 which are gaussian, we obtain, from Eq.( 



X exp 



{pp - g)^ 

I2al 



u. 



D 



Vp exp —iqoT + iq ■ L 



g2 _ ^2 _|_ 

[{Ep - go) - (pp-q) ■ vpf 



l2crlXp 



(10) 



(Pd + q) [{Ed + go) - (pd + q) ■ vd 



i2n 



where 



Xp = ^{v^ + v'^ + vf,]-v 



fp 



'ID 



We now face the problem of performing the integration over g. In usual calculations 
of the processes occurring through the propagation of a virtual intermediate particle the 
integration over its four-momentum g is easily simplified by the Dirac 5-functions arising 
from energy-momentum conservation in the interaction vertices. On the other hand, in 
the wave packet treatment of the initial and final particles, energy-momentum is not 
exactly conserved and there are no Dirac 5-functions available for the simplification of 
the integration over g. However, in our case the production and detection interactions 
are macroscopically separated, so only the propagation of real neutrinos contribute sig- 
nificantly to the process. This physical fact allows us to perform the integration over go 
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by closing the integration path in the lower half of the complex plane. In fact, if the con- 
tribution of the additional path in the lower half of the complex plane can be neglected, 
this procedure picks up only the contribution from the neutrino pole which lies inside the 
path, whereas the antineutrino contribution is neglected. However, the choice of the con- 
tour needs some caution because the term — (go)^ in the exponent diverges as qo ^ —ioo 
and this prevents us from using the usual half circle contour that encircles the lower half 
of the complex plane. Instead, the appropriate integration path has the form of a rectan- 



T. The 



gle whose lower side dissects the imaginary axis at Qq = —i Qa^XpXD / {\p + Xr,) 
contributions from the three sides except the real axis are negligible. The integrals along 
the sides at ±00 vanish since the sides are finite in length and are damped by the —{qoY 
term in the exponent. The lower side gives a finite result with a exponentially damping 
factor — "^apXpXo I {\p + \d) (easily obtained with a saddle point approximation) 
which suppresses strongly its contribution for macroscopic time separations. Therefore 
the integration over go is dominated by the neutrino pole which lies inside the integration 
contour. The resulting amplitude is given by 

where Ea{q) = + tTi^ and 

Q - (lP-jf,[{Ep-Eaiq))~{pp-q)-Vpf 
-^"W = — TWZ^ ^ 



12a2 l2a^Xp 



{PD + q)' [{En + EM) - jPD + Q) ■ VP? ^^^^ 
12a2 l2a^pXD 

Since ap is small and Xp and Xd are of order of unity, the integral over dq is dominated 
by the minimum of Sa{q), which occurs at q — given by 

[Ep - Ea-vp- (pp - qa)] - (pp - qa) + 



A 



p 



+ [Ed + Ea-VD- (pD + qa)] - {pD + qa) = , (14) 



where Ea = Ea{qa) = + t^I and Va = Qa/Ea are the velocities of the mass eigenstate 
neutrinos propagating between the two interaction vertices. A saddle point approxima- 
tion of the integral over dq leads to 



A^J2Ka ^ea Aa CXp -lEa T + IQa ■ L - Sa {%) " ^ " ^"T) Q'' (L - V^T 



a 



where 



A - 1 j^ ^'^Eg-i-qg 

/Det(l]„^ Ea 



(15) 



(16) 



In Eq.([T5|) \/My denotes J2ij ^jMjjV^- for arbitrary vector V and matrix M. 

The amplitude ([15D describes the process under consideration with the assumptions 
that the wave packets of the external particles are sharply peaked around their aver- 
age momenta and the production and detection processes are macroscopically sepa- 
rated in space-time. The amplitude contains the space-time dependent phase factor 
exp —iEa T + iqa ■ L which gives rise to the conventional neutrino oscillations. The ex- 
ponential damping factor exp [~Sa {qa)] implements the overall energy-momentum con- 
servation only within an uncertainty cXp (if \Ep + Ed\ ^ cXp and Ipp + Pnl ^ cTp then 
exp [—Sa {(Ja)] — 1)- Due to the damping factor exp — ^ (^L — Va (^L — Va , 

since the matrices fla are proportional to l/cri ~ cr^ (see Eq.(|l6D), the amplitude in 



Eq.(15) is non- vanishing if the velocities of the mass eigenstate neutrinos satisfy 

L-VaT 



^CT,. (17) 



However, if the mass difference between the mass eigenstates Ua and Ub is such that 
\va — Vb\ T > ax, then the condition (^) cannot be satisfied by both mass eigenstates. 
In this case, at a given time T the amplitude has two (or more) separate peaks in 
space corresponding to the two (or more) mass eigenstate neutrinos and the experiment 
measures only a constant probability (in space) for the flavor changing process under 
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consideration. This is due to the fact that the wave packets of the two mass eigenstates 
are separated by a distance larger than their width and, since they do not overlap, the 
interference term that produces the neutrino oscillations is damped out. 



Cross Section 



The cross section for the process Eq.(0) is given, from Eq.(|T5D, by 

(T(L, T) OC j dP^Aa Al U*^ Uea U*^ CXp [-Sa {%) " (gt)] 
a,b 

X exp 
X exp 



-i {Ea - Eh)T + i -Qb) ■ L 
■ 1 



L-v^t) {l-v,t)--(L- V, t) Q,' (L - V, T 



1^ 



where / dP represents the integration over the 3-momenta and the sum over the spins 
of the final particles; one must also include appropriate average over the 3-momenta and 
spins of the initial particles which are not measured. 

In a practical experimental setting, the distance L is usually a fixed and known 
quantity, whereas the time T is not. Therefore, the cross section at a given distance L 
is given by the time average of a{L,T). We take L along the z direction and integrate 
over time to obtain 

-1/2 



a{L) OC j dPY^AaAl U^a [va^a^Va + Vb^l^V^ 

a,b 

X exp [-Sa (qa) - Sb (qb)] 



X exp < t 




qhz 



^b' 



Va^a ^^a + Vb^b ^^b 




(19) 



1 {Ea - Eb 



2 Va^^ Va + Vb^b Vb, 
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The first exponential term in Eq.(|T^) guarantees energy-momentum conservation 
within the accuracy of the uncertainty principle. The second gives rise to the neutrino 
oscillation in terms of the distance L. The third is a damping factor which describes the 
coherence of the process allowing significant contributions only from the propagation of 
the mass eigenstate neutrinos in the z direction. In fact, since ~ o"^ and L'^/cy'l. ^ 1, 
the integration over the 3-momenta of the final particles receives its dominant contribu- 
tion when this damping exponential factor becomes maximal. For each pair a and 6, this 
maximum occurs when both Va and Vi, are in the z direction. From the stationary equa- 



tion (|Tj), this is realized when all pp, tTp, and vd are in the z direction. In this case 
the matrices Via and fif, are diagonal (see Eq. ([T6|) ). Let us denote with underlines all the 
quantities evaluated at the maximum of the damping exponential and perform a saddle 
point approximation of the integration over the angular variables that parameterize the 
deviation of Va and Vb from the z direction. As shown in Appendix the result can be 
written as 



a{L) (X ^ J dPY,F,,{L)A,SKa^eal^f.bKb[ia^]^'ia + Vi^rb% 
a,b 



X exp 



A/2 






Vaz + 




Ubz 




zz 




zz 




zz^^ + 




.A 



L 



(20) 



where F_^i,{L) is a factor which is weakly dependent on L. Since Kab{L) becomes constant 
for large L, as shown in Appendix we shall neglect F^i,{L) in the following. In Eq.(|20|), 
/ (iP represents the remaining integrations over the momenta of the initial and final 
particles. The factor represents the geometric decrease of the neutrino flux due 

to the distance of propagation, L. It is important to point out here that the second 
exponential in Eq. (pOD cannot be factored out to derive the oscillation probability as in 
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the usual treatment of oscillations. 

Equation (^) also contains a damping factor which decreases exponentially with 
and measures the coherence of the contributions of the different mass eigenstate neutrinos. 
The coherence length for a 7^ 6 is defined by 



T coh 



2,2 



\2 ' 



\ {Vaz - Ubzf {Uaz " UbzY 

beyond which neutrinos do not practically oscillate. This coherence length can be very 
large in the case of relativistic neutrinos, for which \v^^ — v^J -C 1. 

The last factor in Eq. (pOD is due to the time integration and suppresses the interference 
of the contributions coming from the propagation of different mass eigenstate neutrinos 
unless \E_^ — Ej,\ < ap, as it should be from energy conservation in both the production 
and detection interactions. 



4 Relativistic Limit 

As we have emphasized, although our general result given in Eq . (|20|) exhibits characteris- 
tics of neutrino oscillations, the oscillation probability could not, in general, be factored 
out. We now demonstrate that this can be done when intermediate neutrinos are ex- 
tremely relativistic. 

Let us consider a process in which all the intermediate mass eigenstate neutrinos are 
relativistic, i.e. rria -C In this case, the momentum g^^, the energy and the 
velocity v^^^ can be expanded as 

laz=l0z+^-z 

m„ 



Ea = %,+^z + 
i2 



■ ' 2q^^ (22) 

1 _ 



Oz 
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where g^^ is the solution of the stationary equation (|T^ in the z direction for = 
and e^j, ~ "^a/O'oz ^ Iqz given by the solution of the stationary equation (p!^ in the z 



direction to first order in '^a/^'oz' '^'^ lowest order in the relativistic approximation the 
space-dependent part of cr(L) can be factorized as <j{L) = P{L)(Tq, where do is the cross 
section for massless neutrinos. The space-dependent probability P{L) is then given by 



a,b I ^£oz 



X exp 



L 



9 9 \ 2 



(23) 



The first line of Eq.(^) gives the usual oscillation probability for relativistic neutrinos 
(which can be obtained from a quantum mechanical treatment of the neutrino oscilla- 
tions [0 0, D, P]). The second line of Eq.(p3|) contains an exponent which decreases 
quadratically with the distance L and measures the coherence of the contributions due to 
the wave packets of the different mass eigenstate neutrinos. S ince Oo — [^o]zz ~ '-^'xi 
the coherence length becomes 

The length, L^°^, is the coherence length for the neutrino oscillations, i.e. the two mass 
eigenstate neutrinos Va and z/^ contribute coherently to the flavor changing process only 
when L -C L'^ab'-i which case the probability oscillates as a function of the distance L. 



The coherence length given in Eq.(24) is the same as that obtained by physical intuitions 
in Ref . ^ |l^ and from a quantum mechanical wave-packet treatment in Ref . . If the 
distance L is much smaller than the coherence length L^^b^, the damping factor in the 
probability (^) becomes approximately unity and hence one obtains the usual oscillation 
probability. On the other hand, for distances L ^ L^°^ the two mass eigenstate neutrinos 
Va and z/f, contribute incoherently to the flavor changing process. 
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LZ = 2n . (25) 



From Eg. (p3[) , the well-known oscillation wavelength L°^'^ is 
so that Eq.(^) can be written as 

^coh ^ £0.^0.0 _ (26) 

Hence, the maximum number of observable oscillations before the decoherence of the 
wave packets of the different mass eigenstate neutrinos is given by H, ||, |l^ 

rcoh Q 

Nos. = j^-'^. (27) 

Notice that A^^osc is independent of the neutrino mass eigenvalues m^. Since we have as- 
sumed that o"p is much smaller than the energies of the initial and final particles involved, 
o"p < and A^osc > 1- 

Finally, in the last exponential factor of the probability ( pSj) we have retained a space- 
independent damping factor. Since 

22 2 2-1 

m„ — m„ — 1 
eaz - ebz + (28) 

the probability (p3[) for a 7^ 6 does not vanish only when L°^^ > a^. This result is due to 
the time integration: if the neutrino wave packets are larger than the oscillation length, 
the interference terms are washed out. 



5 Conclusions 

In order to illustrate the quantum field theoretical treatment of neutrino oscillations 
without introducing the concept of weak eigenstates, we have discussed a specific flavor 
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changing process (see Eq.(0)) in which the external particles are described by wave pack- 
ets and the mass eigenstate neutrinos propagate between the production and detection 
interactions which are macroscopically localized in space-time. We have calculated the 
time-averaged cross section which is the measurable quantity in the usual experimental 
setting where the distance between the production and detection interactions is known 
but the time separation is not measured. We have pointed out that in general, it is 
not possible to factor out of the cross section a space-dependent oscillation probability 
because the dynamics of the production and detection interactions is not the same for 
the different mass eigenstate neutrinos. However, we have shown that in the extremely 
relativistic limit the usual quantum mechanical oscillation probability can be factored 
out of the cross section. 
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A Wave Packet 



A gaussian wave packet in the momentum space is given by 



27r(Tr 



-3/2 



exp 



ip- X + iE{p)T 



(29) 



where E{p) = y^fp + and Up is the width of the wave packet for simphcity assumed 
to be the same along the three directions. 
In the coordinate space, we have 



^{x,t;X,T,{p)) = J -^^^{p;X,T,{p)) e 



ip-x—iE{p)t 



(30) 



(27r)3/2 

Since the gaussian wave packet in the momentum space is peaked around the average 
momentum {p), neglecting the spreading of the wave packet, one can approximate 



E{p) ~ {E) + v{p-{p)); 



{E) = E{{p)) = ^{p)' + m^ 



(31) 



dE 



Hence the wave packet in the coordinate space is found to be 



i;{x,t;X,T,{p)) 



-3/2 



exp 



i{p)-(x-X)-i{E){t-T) 



(x-X) -v{t-T) 



4(72 



At time t — T the wave packet is peaked at x — X with a width ax given by 

1 



f^x crp = - 



(32) 



(33) 
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in agreement with the uncertainty principle. The wave packet propagates in space with 
group velocity v. 



B Saddle Point Approximation 

The 2-dimensional integration over the angular variables ^ = (^1,^2) that parameterize 
the deviation of Va and Vf, from the z direction in Eg . ([19|) is of the type 

diexp[zA{OL-B{OL'-C{0] ■ (34) 

For large L the integral gets its dominant contribution from the minimum of B[^), which 
occurs for ^ = ^. We expand all the terms around this minimum: 

(35) 



i=l,2 
1 



m) = B + - y: 

i=l,2 



After a change of variable — the integral in Eq.(|34D can be written in 

gaussian form 



exp 



X exp 



iAL - - C 



B 



1 1-1 



(36) 



The final result is 



27r F{L) 
I2 



Det B 



: exp 



iAL - BL^ - C 



(37) 



with 



F{L) = exp 



B 



(38 



From Eq.(R3) it is clear that the space dependence of FiV) is negligible for large L. 
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